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We reconsider in this work the effects of an energy absorption term in the evohition of primordial 
black holes (hereafter PBHs) in the several epochs of the Universe. A critical mass is introduced as 
a boundary between the accreting and evaporating regimes of the PBHs. We show that the growth 
of PBHs is negligible in the Radiation-dominated Era due to scarcity of energy density supply 
from the expanding background, in agreement with a previous analysis by Carr and Hawking, but 
that nevertheless the absorption term is large enough for black holes above the critical mass to 
preclude their evaporation until the universe has cooled sufficiently. The effects of PBH motion 
are also discussed: the Doppler effect may give rise to energy accretion in black-holes with large 
peculiar motions relative to background. We discuss how cosmological constraints are modified by 
the introduction of the critical mass since that PBHs above it do not disturb the CMBR. We show 
that there is a large range of admissible masses for PBHs above the critical mass but well below the 
cosmological horizon. Finally we outline a minimal kinetic formalism, solved in some limiting cases, 
to deal with more complicated cases of PBH populations. 

I. INTRODUCTION A. The critical mass 



The possibility that PBHs may form in the early uni- 
verse is allowed by many theoretical models. Several pro- 
cesses might lead to their formation; for instance, the 
collapse of primordial density fluctuations could have an 
important role in their creation. A detailed review of 
these mechanisms and of their effectiveness for creating 
PBHs can be found in Ref. [1] . Since the density to which 
matter must be compressed in order to create a black 
hole grows with mass, the powerful pressures counteract- 
ing the compression must be taken into account. As a 
result, black holes with mass m << Mq can not form 
in the contemporary universe. Zel'dovich and Novikov 
[2] and Hawking [3] hypothesized that black-holes could 
have been produced at some early stages of the universe. 

PBHs have been studied intensively since the discov- 
ery of Hawking's quantum evaporation, which scales as 
the inverse of the squared PBH mass. Zel'dovich and 
coworkers analyzed the question of the environment ef- 
fect on the individual PBH evolution. Here, we shall re- 
turn to this point and discuss the PBHs behavior taking 
into account the energy content of medium, their proper 
motions and the causal constraints due to existence of a 
cosmological horizon. Since the abundance of these ob- 
jects may be probed experimentally, the consideration of 
the actual evaporation history is in order. We discuss the 
regimes of evaporation and accretion in Section 2. Sec- 
tion 3 presents an asessment of the kinematical effects 
on the PBHs. Section 4 is devoted to a brief discussion 
of a more general formalism to deal with a generic PBH 
population. The conclusions are presented in Section 5. 



II. EVAPORATING AND NON-EVAPORATING 
REGIMES OF PBHS 



The equation that describes the evaporation of black 
holes due to the Hawking radiation is [3] 

dm A(m) . . 

~dt = W ' { ' 

where A(m) counts for relativistic particles emitted 
from a black-hole with mass m. Numerically, A(m) — 
5.3 x 10 25 g 3 s _1 for black holes with masses > 10 17 g and 
A{m) > 7.8 x 10 26 g 3 s _1 for black holes with masses 
< 10 15 g [4]. However, this expression is strictly valid for 
black-holes immersed in a vacuum, since it does not take 
into account the energy content of background. The next 
step is to consider a term which explicitely contains the 
effective energy density and the appropriate cross sec- 
tion for gravitational capture of relativistic particles of 
this background. These effects have been first considered 
by Zel'dovich and Novikov [2] (see also Ref. [5] for a dis- 
cussion) in the past. We reexamine the physical picture 
emerging from the introduction of an absorption term in 
the present work. We wish to take into account relativis- 
tic effects in this description, so that an effective mass 
density can be defined as 

Peff = P + Mp)/c 2 ■ (2) 

If the environment is assumed to be isotropic and ho- 
mogeneous, the cross-section for relativistic particle ab- 
sorption is proportional to the square of the mass 

crym) ~ lo7T — - — . (3J 

c 4 

This expression should be interpreted as an order-of- 
magnitude estimation for the cross-section, as is well- 
known [2] the corrected relativistic expression is in fact 
larger by a factor 1.68 and we shall neglect this small 
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difference in the rest of the paper. However, it can be 
checked the the results depend weakly on the precise nu- 
merical coefficient of a(m) (see for example Ref.[2]). The 
amount of mass-energy absorbed by the black-hole in a 
time interval dt is therefore 



dm _ WnG 2 
~dt ~ \ 



■ ™ Peff ■ 



(4) 



Then, a more complete differential equation for a de- 
scription of the PBHs mass variation is 



dm 
~dt 



A{m) 16ttG 2 



m* 



+ 



-m p eff . 



(5) 



where the background equation of state enters through 
Peff- A critical mass is obtained by setting eq.(5) equal 
to zero and reads 



m c = 



16irG 2 p 



eff 



1/4 



(6) 



A max , the max- 



where we have chosen arbitrarily A 
imum possible value. Any astrophysical or PBH whose 
mass stays above m c at a given cosmological time is able 
to accrete more mass-energy than the amount that evap- 
orates via quantum effects. We will now discuss the be- 
havior of this parameter during the Radiation and the 
Matter-dominated Eras. 



300 GeV) ~ 106.7 associated to the maximum number 
of particle degrees of freedom; and imposed a maximum 
value of A max = A(T B h > 300 GeV) when the black 
hole emits all the known elementary particles. 

A plot of m c for this Era is given in Fig.l. The 
mass contained in the cosmological horizon defined as 
m h = V geo (t)p rad (t), where V geo w \-kt\ is much larger 
than m c except at the Planck epoch. For a model with 
il = 1, rh is given by the usual Th = R(t) J Q cdt' /R(t'). 

Since in the Radiation-dominated Era m,h ~ 4 x 
10 38 (i/l s) g, then sub-horizon PBHs can form and stay 
above the critical mass. By comparing mu and m c we 
find a large range of admissible PBH masses in the ac- 
cretion regime. 

It should be noted that the critical mass expression 
displayed above is strictly valid if the PBHs are scarce 
relative to radiation, or ppbh <C p ra d- If this is not 
the case, the evaporation of these objects may inject en- 
ergy into the medium, which in turn modifies the critical 
mass. A set of equations for a more complete description 
of the latter situation is given by 



H 2 



dm 
~~dt 



8irG 



k 



(Prad + Ppbh) — 



32ttG 2 



A(m) 

H ; — Pradm 



(10) 



[ID 



m'- 



B. Evolution in the Radiation-dominated Era 

In this epoch, the relationship between p ra d and p is 
given by p — p ra dC 2 /3, then the critical mass becomes 



A r , 



m. 



(:'.2: . 



1/4 



(7) 



Numerically, we can express m c in terms of the 
background temperature (in MeV) and the relativistic 
degrees of freedom of the background particles G(T) 
through the expression p ra d = G{T) T 4 yielding 



2.5 x 10 G(T) 



\MeV J 

1/4 



7.8 x lO 26 ^" 1 



9 ; 



(8) 



or, using the well-known time-temperature relation 
T = 1.55 x G(T)- 1 / 4 (t/ls)- 1 /2 Me y 



4 x 10 17 | — 
1 s 



1/2 



7.8 x 10 2 Vs _1 



1/4 



(9) 



and we have assumed the validity of the Standard 
Model, for which one gets a maximum value G max {T = 



^Uprad + iHp pbh + p rad + p pbh 



Q(m) 



(12) 



Where Q(m) is the input power due to evaporating 
black-holes of a given mass, which obviously vanishes for 
m > m c . We consider a flat model of universe without 
cosmological constant. A simple algebraic manipulation 
of these equations yields 



m c oc 



(1 



Ppbh > 
Prad 



1/4 



(13) 



This expression reduces to the former when 
Ppbh/ Prad < 1, as it should. 

In the general case the initial spectrum can be quite 
complicated [6]. In the specific case of scale-invariant 
density perturbations (satisfying the condition 1/3 < 
ZS.pl p < 1) the formation of a continuous PBH mass 
spectrum is expected [7] . If the mass spectrum of PBHs 
favors light objects, a large fraction of the population will 
be evaporating in the Radiation-dominated Era, more- 
over, if ppbh ~ Prad, we may get a significative reheat- 
ing of the background (see, for example Ref.[8]) and the 
critical mass would be diminished in turn. Generally 
speaking the radiation from the evaporating subpopula- 
tion could feed for the non-evaporating massive black- 
holes above m c . Consequently, the spectrum would be 
distorted in a non-trivial way. 
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C. Evolution in the Matter-dominated Era 

Following the same reasoning than in the former point 
it is straightforward to obtain the evolution of critical 
mass in the Matter-dominated Era. We know that the 
scale-factor evolves as R(t) oc i 2 / 3 and the background 
temperature scale as T(t) oc then, using again 

the above formulae, the critical mass for this epoch is 
found to be 



.//,. = 7 x ir'| — | q 



t 



2/3 



(14) 



Where tn denotes the Hubble time of expansion. This 
result obviously applies to astrophysical black-holes as 
well, as all of the latter objects are far above the present 
critical mass and therefore will be accreting photons for 
a long time. One may wonder about the actual amount 
of accretion/mass gain [2,5]. We can check that the rate 
of growth of the critical mass m c is in fact much larger 
than the rate of growth of any particular black-hole from 
a reasonable formed spectrum. In this case, we can set 
w(tff) = mo ~ constant and obtain the approximate 
equality 



mo 



m c (t r , 



(15) 



which is the minimal mass of black holes that survive in 
a closed universe (here, t max = (n/2 H ) ~ 7rf£r/2 
is the cosmological time at maximal expansion of a uni- 
verse with fl > 1). Numerically we obtain too ~ 10 26 g 
for this minimal mass. With this result, we conclude 
that all astrophysical black-holes should survive in closed 
universes. From the maximal expansion, the subsequent 
contraction phase will be raising the temperature of the 
background and these black-holes begin to grow even 
faster than in the previously expanding phase. 



D. Absence of relativistic growth of PBHs 

Now, we will proceed to develop an argument to show 
that similarity solutions do not exist in the Radiation- 
dominated Era. Our arguments are based only on the 
nature of the equation of state for this epoch and the 
existence of a cosmological horizon. By calculating the 
ratio 



= | dm abs /dt 



dm eva p j dt 



rn 
rn r 



(16) 



we imediately conclude that a black-hole with to 3> 
to c absorbs a net positive amount of radiation. Such an 
object may be easily created with a mass between m c and 
nih for all cosmological times (see eq.9). The question is 
: is a relativistic growth of these objects possible? For 
those black- holes with to » m c , the evaporation term 
can be neglected and we get simply 



dm _ 32ttG 2 
~~dt ~ ~ 



■ Prad TO . 



(17) 



On the other hand, in the relativistic growth regime 
the mass increase satisfies 



dm 
lit 



c 

2G 



(18) 



The condition eq.(18) is necessary and sufficient if we 
ignore the transitory existence of a vacuum around the 
black-hole immediately after its formation. Comparing 
eqs.(17) and (18) we get a lower bound to the ambient 
density for a relativistic growth of PBHs to be possible, 
namely 



> 1.2 x 10 



r,2 



TO 



I0 15 g 



g cm 



(19) 



Relating the density to the cosmic time we find that 
the universe can be dense enough to feed a relativistic 
growth of a PBH with a given mass for times earlier than 



trel(m) < 6 X 10 



TO 



(20) 



However, the very existence of a cosmological horizon 
imposes the condition 



t causal > 5 X 10 



-24 ( m 



I s ; 



(21) 



for any causal formation of a PBH of a given mass 
to. Comparing eqs.(20) and (21) we obtain a small 
time interval for considerable accretion to occur, if at 
all. The subsequent expansion quickly diminishes the 
absorption term and the absorbed mass-energy becomes 
a very tiny fraction of the original mass. Then, we con- 
clude that causality considerations leave no room for a 
substantial growth of PBHs due to an increasingly rar- 
efied background. We note that dm/dt — ► c 3 /2G only 
if to — > to/j. Since we expect very massive black- holes 
to be rare in realistic mass-spectra, then, the above ef- 
fects rule out a substantial mass gain for those black- 
holes satisfying m c < to <C to^ and the mass in- 
crease remains a moderate one [5]. Another way to 
understand this effect is to integrate eq.(17) to give 
m(t) ~ too (1 + 3.5 x 10 _39 [TOo/(/][to/s] _1 ); where to is 
the initial mass and to the formation time of the black- 
hole in the Radiation Era. Since the gain for black-holes 
inside m c < too <C mu is very small, we again conclude 
that a relativistic growth is precluded. 

It is important to notice that, since the critical mass 
itself grows much faster than the mass of those black- 
holes that stay above it, we can estimate the time a 
given black- hole (formed at tF with original mass too = 
to(£_f), above m c by hypothesis) spent in the accret- 
ing regime inside the Radiation-dominated Era. This 
period is denoted by t c (m ) and is given by condition 
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m(t c ) — m c (t c ). Since r = (m/m c ) 4 , we may ignore the 
Hawking term when mo 3> m c (tp), and invert an ana- 
lytical solution of the condition m(t c ) = m c (t c ), yielding 



tc 



400 



Vio 15 s 



(22) 



in the Radiation-dominated Era and t c ~ 2 x 
10 10 tn (m/m ) 3 / 2 in the Matter-dominated Era. We 
may compare this value to the evaporation time due 
to quantum effects alone tuawking — 10 17 s(m/l 15 g) 3 . 
We conclude that the evaporation time picks up a small 
correction due to the background influence, which delays 
the onset of the evaporation regime unless the black hole 
is more massive than 1 23 g. In the latter case the non- 
evaporating regime has a duration comparable to the age 
of the universe. 

In spite that these considerations show that a PBH 
with, say, m > 3 x 10 22 g slowly stays accreting en- 
ergy during the whole Radiation-Dominated Era, the net 
gain of energy is very small. For example a PBH formed 
at t ~ Is with m ~ 10 20 g gains a mass fraction 
of Sm/m 10~ 19 when accreting for 10 s s. We have 
checked that this feature holds irrespective of the actual 
relationship between tp and m . Then, all those black- 
holes with m >> m c are better described by dm/dt = 
provided they form a rarified gas. If the transient vac- 
uum formed around the black-hole immediately after its 
birth is considered, accretion is impossible before a time 
of the order of H , and the background radiation would 
be even more rarified than our previous estimates when 
the steady conditions are re-established. 



III. KINEMATIC AL EFFECTS 

Up to now we implicitely assumed that the PBHs are 
at rest relative to the thermal background or have negli- 
gible peculiar velocities, However, the primordial events 
ocurring in the early universe might be responsible for 
imparting substantial kinetic energy to black-holes. For 
instance, if a randomly moving black-hole is hit by an 
expanding bubble wall, it is possible that the former can 
acquire a large momentum because of the sudden depo- 
sition of energy from the relativistically moving bubble 
wall [6] . A black-hole moving with velocity v pec sees the 
background temperature to have an angular dependence 



T 



7 (1 



■cos 0) ' 



(23) 



with T being the temperature in the cosmic frame. 
Due to this peculiar motion, the rate of change of mass 
of the black-hole in its rest frame is given by 



dm 
lb 7 

where 



A(m) 32ttG 2 



< Prad > m 



(24) 



<prad>= J dflp rad (T(^,e)) 



(4 7 2 - 1) 



Prad •> 



(25) 



is the average over the angles and r is the proper time. 
After a simple calculation we get, using the relation be- 
tween the proper time and the cosmic time 7 dr = dt 



dm _ A(m) + 32ttG 2 /4 7 2 - 1 
dt 7m 2 c 3 \ 37 



Pradm 2 . (26) 



A high value of the Lorentz factor 7 enhances the 
energy-mass absorption term, and at same time inhibits 
the mass loss due to the Hawking radiation. These ef- 
fects may be incorporated in the critical mass definition, 
which becomes 



m c (T, 7 ) = 



4 7 2 



-1/4 



m c (T) ; 



(27) 



where m c (T) is defined in eq.(8). 

How do we interpret these results? For an observer 
moving with the black-hole, the thermal background is 
not isotropic and the mass accretion can be larger than 
the evaporation. The critical mass is decreased due to 
the Doppler effect. From the point of view of an asymp- 
totic observer, comoving with the cosmological expan- 
sion, both terms are altered: the Hawking radiation term 
is reduced by the Lorentz factor due to the time dilation 
effect (in the rest frame the Hawking radiation depends 
on the rest mass only). The second term is altered due 
to same Doppler effect and also corrected by the Lorentz 
factor. If we have to decide whether a particular black- 
hole is above or below the critical mass, taking into ac- 
count its possible peculiar motion, we must resort to the 
rest frame or the comoving frame. Any other frame is 
pcrtubed by the proper motion of the observer and the 
results are altered. For a black-hole with rest mass m 
we can calculate the minimum peculiar motion necessary 
to ensure that the absorption term overcomes the Hawk- 
ing radiation term. Thus, we can define the threshold 
Lorentz factor 



f 



(28) 



A black-hole which is above this bound (see Fig. 2) will 
accrete mass, even if m < m c , and may stay in this con- 
dition for a long time depending on the actual evolution 



of the v pec field. 



IV. AN OUTLINE OF A KINETIC THEORY 

In order to study the detailed evolution of a PBH pop- 
ulation one must go beyond the thermodynamical de- 
scription and set up a formalism capable of describing 
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the variety of physical effects existing in the early uni- 
verse. Such an analysis may be achieved through a set 
of differential equations for the mass distribution func- 
tion of black-holes /, the cosmological expansion, and 
the equation for mass variation. The initial conditions 
must also be specified (for example, the fluctuation pat- 
tern that originates a given initial mass spectrum). A 
minimal kinetic formalism is given by 



dt 



dpdf_ 
dt dp 



dt 2 dfi 



HB df 

Hp op 







H 



8nG k 

g \Prad + Ppbh) ~ ~j^2 



dp 

~dt 



A(m) 
m pl [i 2 



32 7T G 2 m P i 2 

H ^3 PradP 



(29) 



(30) 



(31) 



where Ro is the scale factor at the formation time for 
this population (assumed to be unique for simplicity). 
Furthermore, since we have assumed a dilute PBH gas, 
a separation of variables is possible to cast f± into the 
form 



f±(p,M 



F±(p,t)J±(p) 



(37) 



Within our hypothesis p p bh <C p ra d (and also that 
the bulk of the PBH population residing between to* = 
(3A/ H) 1 / 3 and m c in order to avoid a substantial injec- 
tion of energy due to low-mass objects), the third term 
in eq.(28) can be neglected and we get 



F± 



dt 



3H 





/8F±\ 


Hp 


fdJ±\ 


F± 1 


V dp J 


J± ' 


\ d/3 J 



0. 



(38) 



4 H p r ad + 3 H ppbh + prad + Ppbh 



Q(m) 

i? 3 C 2 



(32) 



where p is the dimensionless PBH mass scaled to the 
Planck mass mpi, / is the distribution function defined 
on an extended phase-space including the mass. With 
these definitions we have p p bh(iTi,t) = mn(m,t). Given 
some initial conditions the solution of this set of equa- 
tions describes a dilute black-hole gas in a non-rclativistic 
regime for a given model of the universe. Now, we will 
seek particular solutions for a flat model of universe, a 
dilute and sub-dominant PBH-gas, and non-relativistic 
peculiar motions. Since that the (^)(§^) term changes 
sign as a PBH crosses the critical mass, we split / in two 
functions defined by 



and 



/+ = f(P > Pc,t) 



f- = f(P < Pc,t) 



(33) 



(34) 



Since it can be shown that there is no termodynamical 
equilibrium between the PBH gas and the background 
radiation, these objects cannot stay fixed at the crit- 
ical mass (mainly because of the cosmological expan- 
sion). The relevant observables can be calculated if f± 
are known. When expressed in terms of f± the comoving 
number density of PBHs is simply given by 



n± 



g{p) 

(2tt) 3 



d 3 P f±(p 1 P 1 t) 



(35) 



with g(/j,) is the statistical weight of the black hole gas 
and the integration runs over all the peculiar momen- 
tum p = mpicpLp. The geometrical expansion of the 
universe can be isolated with the choice 



/±(m,/M) 



ffi 



f±(p,M 



(36) 



The quantity dj3 1 J± (/?) is the comoving number of 
black holes with peculiar velocities belonging to the in- 
terval (/3, (3 + d/3). We will show below that the minimal 
kinetic formalism is consistent with the concepts deduced 
in the previous sections, in particular the negligible mass 
accretion for black holes described by /+. 

The number density of the evaporating n_ and non- 
evaporating n+ populations is given in terms of the new 
variables as 



n±(p,t) 



g(p)(m P icf 
2tt 2 



P 3 F ± 



dp f3 2 J±((i) ; (39) 



respectively. The kinematical effects discussed in sec- 
tion III can be taken into account by replacing the up- 
per limit of the integration over J_ by Pth(p) = — 
l/lth), where jth is given by eq.(27). We shall neglect 
this effect in the following. 

We finally arrive to the expressions of the comoving 
mass density for each subpopulation as 



p(p > p c ,t) = 
and 

p(p < p c ,t) 



_ g{p) m pl c 3 1+ 



2tt 2 



J fi 



Mh(*F) 



(*) 



g(p) '• 



"pi • 



2tt 2 



/■jUc(t) 

Jo 



dpp 3 F+(n,t) , 
(40) 

/'/' '/•' •:/'•':■: 
(41) 



where we have extended the upper limit of the first in- 
tegration since there is no mass gain for black holes with 
initial mass too < mh(tp), as discussed in section II. D ; 
and defined I± = ^ dp p 2 J±{P). 

A generic PBH population will show one novel feature: 
if there are some black holes with too > m c (tp) these 
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objects will be crossing the critical mass curve and go 
into the subcritical region. The kinetic formalism is able 
to take into account this effect in the following way; let 
us separate p(p < p c ) as a sum of two contributions p\ 
and p 2 , defined by 



pi(p < p c ) oc 



dp p 3 F_((j,,t) 



(42) 



and 



Mc(t) 



P2(p < Pc) oc / dpp 3 F_(p,t) (43) 

where p*{t) ~ p c (t F ){l - 3A{t - t F )/mVc(tF) 3 ) 1/3 
is a boundary curve, which evidently divides the evap- 
orating black-holes into a population that was already 
born in such condition, and a population that entered 
the evaporation regime after crossing m c in a finite time 
(see Fig. 3). Since the latter PBHs satisfy 5m/m <C 1, 
then their number density is described by the usual di- 
lution of non-relativistic objects n = n (R/R )~ 3 and 
therefore we get 



rm h (t F ) 

p(p > p c ,t) = / dmn(m,t) 

Jm c (t) 



= n 



R(t) 
Ro 



(m h {t F ) - m c (t)) 



(44) 



This density is diluted by two factors, as expected (see 
Fig. 3): the first is the usual geometrical expansion and 
the second their migration into the subcritical regime. It 
is obvious that p(p > p c ,t > t c (ph)) = 0, since after 
tc(fih) no PBHs in the accretion regime remain. Com- 
paring eqs.(39) and (43) we get that F + cx p~ 3 . 
Substituting F + into eq.(37) and solving for J + we ob- 
tain 



j+(I3) = j /r G+(t) 



(45) 



where the exponent is G+(t) = 3(1 + 2p/pH) (with 
p/pH <C 1) and Jo is defined by the initial conditions. 

On the other hand, the subcritical population is much 
more complicated than the supercritical one. To sim- 
plify the solutions we assume that the source term in 
the last eq.(31) may be written as Q(m)/R 3 (t)c 2 - 



p(m) / T evap (m) , where 



evap 



(m,Q — m 3 )/3 A(m) is 



the timescale for a black-hole to evaporate from mo to 
m, and get the following solution 



m-3 

F- oc t-j— (exp(-tri) - p*) + 

+ ex P (-tr 2 )e(Ai» - p)) 



(46) 



where is the step function and l/T\, l/r 2 are mass- 
dependent timescales defined by 



3A(m) 



Pa < Pc(t F ) 



(47) 



and 



T, = 



(j4 - H 3 )m 3 Pl , 

-5—77 — r— 7 — r- Pa > Pc(tF) 
3A(m)t c (p ) 



(48) 

Therefore, substituting F_ in eq.(37) we obtain the 
solution 



J-ifi) = J /3- G - (t) ; 
and the exponent is given by 



G-(t) 



3 1 



X 



dX 



2p 

~p~W ~ 3H X ~ 3XH~dp~ 



(49) 



(50) 



where X, etc. are time derivatives and X is defined as 
X = cxp(-tTi) Q(p - p*) + cxp(-tr 2 ) Q(p* - p) . 

(51) 

The exponential factors exp(— tYi) explicitly describe 
the evaporation of PBHs. In fact, Ti controls the rate of 
evaporation of PBHs from p to p, and T 2 is the anal- 
ogous timescale, but corrected for the time t c (po) when 
no evaporation occurs (see Fig. 3). 

An important point (not addressed in this work) is to 
relate the mass-spectrum functions F±(p,t) to the ini- 
tial conditions which formed the PBH population, for 
example, large primordial fluctuations imprinted in the 
radiation field. A handful of works in the literature ad- 
dress this problem, i.e. Ref.[8]. Carr's mass spectrum 
corresponds to the particular choice F = p n ~ 3 N(p, t) 
(where N(p,t) = Ap~ n and A — constant) in the pre- 
sented formalism. The accreting subpopulation of Carr's 
spectrum F + is not specially affected by our considera- 
tions because that each of that objects remain near the 
formation mass. However, the evaporating part F_ de- 
pends quite strongly on time, and therefore the initial 
spectrum is modified with yet unknown consequences. 
Incidentally, eqs.(45) and (49), show that the the PBH 
population does not enter into equilibrium with the am- 
bient radiation in this approximation. 

As an application of the above results we may estimate 
the maximum allowed number density (at a given forma- 
tion time) consistent with the observed expansion and the 
present bound Q p bh < 10~ 8 [9]. If we consider a power- 
law mass spectrum of the type N oc p~ n and a "burst" 
of PBHs formation at t = tp'm the Ratiation-dominatcd 
Era, the above expressions can be manipulated to give 

n(t F ) < 3.5 x 10- 42 ^ n_1) (^) 



t F 



-5/2 



x| — j (1 + 6/2) n pbh h 1 cm~ 



(52) 



with S = m c (t h )/m h (t F ) < 1. 

This discussion illustrates some features of the for- 
malism, which is of course not restricted to Radiation- 
dominated, flat cosmological models. More general con- 
ditions and a variety of mass spectra will be presented in 
future works. 
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V. CONCLUSIONS 
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VI. FIGURE CAPTIONS 



Figure 1. Critical mass m c as a function of the temperature in the Radiation-dominated Era. The curve splits 
this plane in two different regimes as throughly explained in the text. Note the jumps due to the changing degrees of 
freedom of the background. 

Figure 2. Threshold Lorentz factor j t h as a fuction of the mass ratio m c /m. Black holes above the curve accrete 
matter and do not imediately evaporate (in spite of the condition m c /m > 1) due to the 7 dependence of both the 
accreting and evaporating terms of eq.(26). Those black-holes satisfying m c /m < 1 do not evaporate even if their 
peculiar velocities vanish. 

Figure 3. Qualitative graphical representation of the evolution of a generic PBH population in the \i — t plane. The 
crossing-time t c (/xo) for a black hole with initial mass no is depicted, together with the horizon mass \Xh and critical 
mass fi c defining the different regimes in the Radiation-dominated Era (see the text for details). 
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